DRIFT OF PHASE FLUCTUATIONS 
IN THE ABC MODEL 
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Abstract. In a recent work, Bodineau and Derrida analyzed the phase fluc- 
tuations in the ABC model. In particular, they computed the asymptotic 
variance and, on the basis of numerical simulations, they conjectured the pres- 
ence of a drift, which they guessed to be an antisymmetric function of the 
three densities. By assuming the validity of the fluctuating hydrodynamic ap- 
proximation, we prove the presence of such a drift, providing an analytical 
expression for it. This expression is then shown to be an antisymmetric func- 
tion of the three densities. The antisymmetry of the drift can also be inferred 
from a symmetry property of the underlying microscopic dynamics. 



1. Introduction 

The ABC model, introduced by Evans et al. [121 HE], is a one-dimensional stochas- 
tic conservative dynamics with local jump rates, whose invariant measure undergoes 
a phase transition. It is a system consisting of three species of particles, traditionally 
labeled A, B, and C, on a discrete ring with L sites. The system evolves by nearest 
neighbor particles exchanges with the following rates: AB —> BA, BC — > CB, 
CA -> AC with rate q and BA -> AB, CB -> BC, AC -> CA with rate l/q. 
In particular, the total number of particles N a , of each species a £ {A, B, C}, are 
conserved and satisfy Na + Nb+Nc = L. When q^l, Evans et al. [12j[l3] argued 
that in the thermodynamic limit L — > oo with N a /L — > r a the system segregates 
into pure A, B, and C regions, with translationally invariant distribution of the 
phase boundaries. In the equal densities case Na = Nb = Nc = L/3 the dynamics 
is reversible and its invariant measure can be explicitly computed. As shown in 
[14t ITS] , the ABC model can be reformulated in terms of a dynamic of random 
walks on the triangular lattice. 

As discussed by Clincy et al. [10], the natural scaling to investigate the asymp- 
totic behavior of the ABC model is the weakly asymmetric regime q — exp { — } , 
where the parameter /3 plays the role of an inverse temperature. With this choice, 
the reversible measure of the equal densities case ta = tb = tq = 1/3 becomes a 
canonical Gibbs measure with a mean field Hamiltonian, which undergoes a sec- 
ond order phase transition at j3 c = 27rv3- More precisely, for (3 < f3 c the typical 
densities profiles are homogeneous while for (3 > j3 c the three species segregate. 
For unequal densities the invariant measure of the ABC dynamics on a ring is not 
reversible and cannot be computed explicitly. The asymptotic of the two-point cor- 
relation functions in the homogeneous phase is obtained in (6] [10] , where the large 
deviation rate function for the stationary measure is also calculated up to order /3 2 . 
When the ABC dynamics is considered on an open interval with zero flux condition 
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at the endpoints, the corresponding invariant measure is reversible for all values of 
the densities pQ. In particular, it has the same Gibbs form as the one in the ring for 
the equal density case. In this paper we shall however stick to the case of periodic 
boundary conditions. 

In the diffusive scaling limit, the hydrodynamic behavior of the empirical den- 
sities p = p(x,t) — [pA{x,t), ps{x,t), pc(x,t)) , where t > and x £ T, the 
one-dimensional torus of length one, is given by the parabolic system, 

d t p a = dxxPA +/3d x [pa(pb - Pc)] , 
< d t p B = d xx p B + /3d x [p B {pc ~ Pa)], (1.1) 
dtpc = d xx p c + Pd x [pc(pa - Pb)] ■ 

The initial conditions arc determined by the starting microscopic configuration, 
they satisfy the constraints p^(-,0) + ps(-,0) + pc(-,0) = 1 and < p Q (-,0) < 1, 
a £ {A,B,C}, which are preserved by the above flow. Clearly, it also preserves 
the total mass of each species, i.e., r a — Jdx p a (x, t), a € {^4, B, C}, is constant in 
time. We mention that in the equal densities case, r a = rg = vc = 1/3, the flow 
defined by is a (suitable) gradient flow of the large deviations rate function, 
which in this case reads, 

Fp{p)= \ dx y2Pa(x)logp a (x) +p / dx dy ^2p a {x)p a+2 {y), (1.2) 
Jo a Jo Jx 

where the sum on the species label a is modulo three. This is the macroscopic 
counterpart of the reversibility of the underlying microscopic dynamics. 

In the equal density case the stationary solutions to correspond to the crit- 
ical points of (jl.2p . which have been analyzed in pQ. In particular, the homogenous 
profile r = (3,3,3) is the unique stationary solutions for /? S [0,/3 c ], while for 
/3 > (3 C there exists a one-periodic inhomogeneus stationary solution, unique up to 
translations, which minimizes the large deviation functional and it is therefore sta- 
ble for the flow In the general case, the homogeneous profile r = (r^, rs,rc) 
is clearly still a stationary solution to (|l.ip . For j3 small enough (depending on 
7') it is the unique one. As discussed in [6j 110) , a linear stability analysis shows 
that for P > (3 r — 2ixj\J\ — 2r 2 , r 2 = r\ + r 2 B + r 2 c , it becomes unstable. As 
stated there, the phase transition, at least for particular values of r, is expected 
to become of first order. Namely, for some /3 £ (0,/? r ) there exist other stationary 
solutions to (jl.l|) which actually describe the typical profiles with respect to the 
invariant measure of the underlying microscopic dynamics. In the stationary 
solutions to are characterized and their explicit expression is written in terms 
of elliptic functions. This analysis reveals that for some r and j3 £ (0,/3 r ) there 
exist one-periodic stationary solutions to in agreement with the conjectured 
occurrence of a first order transition. 

The next natural issue on (II. lj) is whether it admits traveling wave solutions, this 
is indeed a side question both in [5] and [11]. A simple argument, discussed in [20] 
shows that such solutions do not occur. For the reader's convenience we report the 
argument in Appendix [X] We emphasize that the periodicity of the space variable 
is crucial. Indeed, when the system (|1.1[) is considered on the whole line it does 
admit traveling waves. 

The main result of the present paper is the identification of the drift for the phase 
fluctuations. In order to describe it, fix values of r and (3 in the low temperature part 
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of the phase diagram. The analysis in implies the existence of a one-periodic 
profile p = (paj Pbi Pc) sucn that the one-periodic stationary solutions to (II. ip arc 
given by the translations of p. Observe that, as discussed in [HE], for larger values 
of (3 there exist also ^-periodic solutions, ^-periodic solutions,... However, as proven 
in [I] in the case of equal densities and suggested by numerical evidences in general, 
these less segregated profiles are not expected to describe the typical behavior of 
the microscopic evolution. Consider now the microscopic dynamics in which the 
initial distribution of the species is associated to the macroscopic profile p. The 
hydrodynamical description discussed above implies that on the microscopic time 
scale 0(L 2 ) the associated macroscopic profile does not move. Since macroscopic 
fluctuations are 0(L~5 ), by taking into account the translation invariance, at times 
0(L 3 ) the macroscopic profile is expected to perform a random motion on the set 
■ — z), z £ T} . This motion is refereed to as phase fluctuations. In agreement 
with the above dynamical picture, we observe that in the equal densities case, by 
sampling the particles according to the Gibbsian invariant measure, the law of 
large numbers for the empirical densities of the three species is p(- — £), where £ is 
a uniform random variable on T, see [2] for a formal proof of this statement. 

Bodineau and Derrida [5] have computed the variance of the phase fluctuations 
by using both the methods of the fluctuating hydrodynamics and of the macro- 
scopic fluctuation theory. Moreover, on the basis of numerical evidences, further 
confirmed in [TT] , they conjectured the presence of a drift which they guessed to be 
an antisymmetric function of the three densities. Observe that for equal densities, 
reversibility readily implies that the drift vanishes. In this paper, by assuming the 
validity of the fluctuating hydrodynamic approximation for the ABC model, we 
prove in general the presence of such a drift, provide an analytical expression for 
it, and show its expected antisymmetry property. We also show that this anti- 
symmetry can be inferred from a general symmetry of the underlying microscopic 
dynamics. 

In order to obtain the statistics of the phase fluctuations, we need to obtain an 
effective dynamics on the manifold {p( — z), z € T} . The fluctuating hydrodynamic 
approximation allows to describe the behavior of the microscopic ABC dynamics by 
the hydrodynamic equations perturbed with an additive noise of order . 
In particular, the behavior of the ABC model on time scale 0(L 3 ) is captured by the 
asymptotics of the fluctuating hydrodynamics on the time scale 0(L). Under this 
assumption, the variance of the phase fluctuations can be computed by considering 
the projection along the manifold of the noise: this is the sum of independent, mean 
zero, order random variables, that converges to a Brownian motion. On the 
other hand, the origin of the drift is much less evident and due to the nonlinearity of 
(jl.ip . Indeed, we show that in a time step T with 1 <r<L the nonlinear term 
gives a deterministic contribution to the effective dynamics of order T/L, which 
sum up to a finite drift at times O(L). 

We finally mention that the statistics of dynamical phase fluctuations has been 
already analyzed in the context of the one-dimensional nonconservative stochastic 
Ginzburg-Landau equation 8 , 9, 16 . In particular, if the reaction term is not sym- 
metric, the mechanism outlined above gives rise to a constant drift in the resulting 
random motion. This has been rigorously proved in [7] . Phase fluctuations for the 
Kuramoto model, which is a mean field conservative dynamics, have been recently 
discussed in [3]. 
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2. Identification of the drift 

In this section we introduce the fluctuating hydrodynamic assumption and iden- 
tify the drift of the phase fluctuations in terms of the semigroup generated by the 
linearization of (jl.ip around the periodic profiles p. We shall frequently refer to [5] 
and use the same notation introduced there. For convenience we set e = l/L. 

The fluctuating hydrodynamic assumption. At the macroscopic level, the 
effect of the microscopic fluctuation can be modeled by adding to the hydrodynamic 
system (jl.ljl a suitable random force whose statistics can be inferred by an informal 
computation on the underlying Markov dynamics. Referring to either Appendix IBl 
or OEJUU] for the details of such computation, the corresponding stochastic system 
reads, using vector notation and denoting partial derivatives with subscripts, 

/ Pa(pb - pc) \ _ ( Va \ 
+ f3 Pb(pc-Pa) +V^\ Vb (2-1) 

V P C {PA ~ PB) ) x V J x 

where t > 0, x e T = R/Z, the one-dimensional torus, and, conditionally on the 
value p, the noise rf — (j]\ 1 'r}%,'rfc) ^ s Gaussian with correlations 

(»£(*> 0) = S a ,a'(P5 x, x') 5 e (x-x') 5(t-t% a, a' e {A, B, C}. (2.2) 

Above, 6 £ is the e-approximation to the Dirac's S- function, i.e., its length scale is 
of order e, and S is the matrix 

E(p; x, x') = D{p; x, x') + D(p; x',x) (2.3) 

in which D(p\ x, x') is 

f Pa{x)\pb{x') + pc(x')] -p A (x)p B (x') -p A (x)p c (x') 

-p A {x)p B {x') p B {x)[p A {x') + p c {x')] -p B {x)p c (x') 

-Pa{x)pc(x') -p B (x)p c (x') pc{x)[pa{x') + p B (x')} 

Observe that E has a vanishing eingenvalue with corresponding eigenspace spanned 
by the vector (1, 1, 1) T . Accordingly, (|2.ip preserves the constraint J2 a P a = 1- 

Clearly, the (deterministic) hydrodynamic system (11.11) is recovered from (12. ip 
simply by setting e = 0. As discussed e.g., in [3TJ § II.3.5], (|2.ip also predicts 
the nonequilibrium Gaussian fluctuation which can be inferred by linearizing (I2.ip 
around an hydrodynamic solution. As discussed in the Introduction, in our analysis 
we shall need the stronger assumption that (|2.ip correctly encodes the behavior of 
the ABC dynamics on time scales longer then the hydrodynamical one. More 
precisely, as phase fluctuations becomes observable (i.e., macroscopically of order 
one) at microscopic times of order L 3 , the connection between the analysis on the 
stochastic system (|2.1I) performed in the sequel and the microscopic dynamics relies 
on the hypotheses that the ABC dynamics at times L 3 can be captured by looking 
at (|2.1I) on times of order e _1 and then taking the limit e — » 0. 

In regard to the fluctuating hydrodynamic assumption, we simply observe that 
the use of nonlinear stochastic equation to describe the evolution of particle systems 
beyond the hydrodynamic scale is a common practice in nonequilibrium statistical 
physics. On the other hand, the rigorous justification of such procedure is a most 
challenging task of mathematical physics, see [3] for the weakly asymmetric exclu- 
sion process which has, albeit much simpler, similar feature to the ABC model. 
We finally emphasize that in (|2.ip we have somehow kept track of the underlying 
discrete structure by using a colored noise with spatial correlation length of order 
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Denoting by 



s. With this choice, we avoid in particular the difficult problem of giving a precise 
mathematical meaning to (|2.I|) with a space-time white noise, see [17] for the KPZ 
equation whose nonlinearity has a similar structure. As it will clear in the com- 
putation of the drift, the problem of ultraviolet singularities will however appear 
when taking the limit e — > 0. 

The effective phase dynamics. It is convenient to take advantage of the con- 
straint pa + Pb + Pc = 1 an d write the system (|2.1[) in terms of 

the nonlinear transport and the noise, respectively, we thus rewrite (|2.1[) as 

+ /3AA(p) ;c + % /i77^. (2.5) 

Conditionally on the value p, the noise is Gaussian with correlations as in (|2.2I) , 
where x, x') is the 2x2 matrix as in (12.31) in which now 

D(p;x,x')= ( PA{X)[ ) - p *W ] 7& )M i\y (2-6) 

Hereafter, we fix the total densities r — (r^rg) and ft > in the low tempera- 
ture part of the phase diagram established in [IT] . The analysis therein implies the 
existence of a one-periodic profile 



P = P(x) 



K PB(X); 

which solves the following system of ordinary differential equation, 

p" + ftAf{p)' = 0, (2.7) 

and satisfies the constraints p a > 0, Pa+Pb ^ 1j an< i /dx p a = r a , a G {A, B}. The 
one-periodic stationary solutions to (|2.5[) with e = are given by the translations 
of p. 

We denote by p z {x) — p(x — z) the translation of p by z on the torus. Let also 
C z be the linear part of (|2.5|) around p z , i.e., 

C z ip = ip xx + ft(B z ip) x , (2.8) 

where 

B -W "»(*-«). B =( % -S _1 -2,,?& B + l)- <"> 
Recall that 

X*0) = (p B (x-2),-p j4 (x-z)), p' 2 (x) = 

are the left and right eigenvectors of C z with zero eigenvalue, which can be easily 
shown to be a simple eigenvalue. The rest of the spectrum is composed by a 
countable set of eigenvalues. As in [5], we assume they have strictly negative real 
part, bounded away from zero. In regard to this assumption, we remark that in 
the equal density case C z can be realized as a nonnegative self-adjoint operator 
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on a suitable Hilbert space and the above assumption is fulfilled. A standard 
perturbation argument implies that it also holds for r close to (i, |, |). Under this 
assumption, the manifold Ai — {p z : z € T} is locally exponentially attractive for 
the deterministic flow (jl.ll) . The (not orthogonal) projection onto the null space of 
L z is given by the tensor product jp' z (x)x z (x'), where 

1 f 1 

-= / dx X P- (2.10) 
7 Jo 

We shall also denote by V z the projector whose integral kernel is given by V z (x,x') — 
S(x — x')l — r )p' z {x)xz{x'). Observe also that the right eigenvector p' z is the infini- 
tesimal generator of the translations on Ai . 

We consider the random flow (|2.1[) with a initial condition lying in a £5 -neighbor- 
hood of Ai. As fluctuations orthogonal to Ai are exponentially damped by the 
deterministic part of the flow, we deduce that the solution remains in such neigh- 
borhood until a large fluctuation takes place. For our purposes, those large fluc- 
tuations can be neglected, as they occur on time scales exponentially large in s . 
In particular, with very large probability, such solution to (|2.1j) remains in a £2- 
neighborhood of Ai up to time scales of order e , which is the relevant one for 
the phase fluctuations. 

In order to describe the motion along the manifold Ai we use the Fermi coor- 
dinates ((,ip), defined as it follows. Given p in a neighborhood of Ai, the angular 
coordinate C € T, called the center of /?, is defined as the point z such that the 
component of p — p z in the null space of C z vanishes, namely, the solution to 

F(z)= f dx Xz (p-p z ) = (2.11) 
Jo 

and then we set ip — p — p^ , so that if) = V(ijj. 

According to the scheme introduced in [9] and further developed in [8], the 
motion along the manifold Ai can be identified by the following recursive procedure. 
Given an initial datum p(-, 0) in a £2 -neighborhood of Ai, we let Co be its center and 
decompose p(-,0) = /%(•) + VK'jO)- We then linearize the evolution (12.51) around 
P(- and compute the displacement in a time interval 1 < T< s . At this point, 
we recenter the solution by computing the new center £r at time T and then iterate. 
As it will be clearer in the sequel, after (eT) _1 steps we get a finite displacement 
of the center, which has the form of a Brownian motion with a constant drift on T. 
Dropping the subscript Co from the notation, we now detail the first step of such 
procedure. By introducing u = u(x, t) as 

( u a\ 
U =[u B )=P-P> 

the evolution (12.51) can be recast into the form, 

u t = Cu + pM 1 (u) x + y/eff x + --- , (2.12) 

where 

C is as in (|2.8j) - (|2.9j) . and we approximated rf with the Gaussian noise fy e , having 
covariance matrix £(/?; x, x') — D(p; x, x') + D(p; x' , x), with D(p; x, x') as in (|2.6p . 
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By Duhamel formula, 

u(t) = e tc V(0) + /3 [ dse^- s ^ c Mi(u(s)) x + ^W e {t) + • • • , (2.14) 
Jo 



where W £ (t) = 




is Gaussian with covariance 



(W e (x,t) W £ (x',t) T ) 

t A rl _ (2.15) 

ds dy dy Gy(x,y,s)5 e (y-y)'E(p;y,y)G y (x ,y ,s) 
Jo Jo 

and G(x,y,t) = e tc (x,y) is the fundamental solution associated to C. 

Since ip(0) has vanishing projection on the null space of £, i.e., 'PV'(O) = -0(0) , 
the first term on the right hand side of (|2.14j) is exponentially small as t — > oo. 
In particular, u(T) is of order e' so that we can compute the displacement of the 
center in the time interval [0, T] by solving (|2.14l) in the linear approximation. We 
thus get 

Ct = Co-7 / dx X (x)u(x,T) + ■ ■ ■ , (2.16) 
Jo 

where 7 is defined in (|2.10l) . 

Computation of the variance and the drift. The variance of the infinitesimal 
displacement Ct — Co can be effectively computed by considering the component of 
the noise W £ in the direction spanned by the right eigenvector p'. Recalling (|2.15|) 
and that x e ' £ = 

((Ct - Co) 2 ) « e 7 2 ( ( J 'dx X (x)W s (x, T))') w e T 7 2 J dy x'ivMP; V, y)x'{y) T 
= ZeT^ 2 Jdy [p A (l - p A ){p' B ) 2 + 2p A p B p' A p' B + p B (l - p B ){p' A ) 2 ] , 

where we evaluated (|2.15p with S e replaced by the true Dirac's (5-function. By 
translation invariance, each step of the iterations has the same variance. Therefore, 
since the fluctuation of the infinitesimal displacements relative to different steps are 
also independent and mean zero, they sum up to a Brownian motion with variance 
a 2 = a 2 ((3,r) given by 

^ = V Jdy [p A (l - p A ){p B f + 2p A p B p' A p' B + p B (l - p B )(p' A ) 2 ], 

which has been computed in [5j Eq. (7)]. 

As discussed in the Introduction, the non-linear term contributes to the phase 
fluctuations at times e _1 , by giving a constant drift to the resulting random motion. 
Our next aim is to identify this drift as a function of (3 and the total densities r A ,r B . 

As for the variance, we discuss in detail the contribution picked up in the first step 
of the iteration. Indeed, again by translation invariance, each step of the iterations 
gives the same contribution to the drift. Since u is of order and the non-linear 
term Mi in (|2.13[) is homogenous of degree two, the latter gives a contribution of 
order eT to the infinitesimal displacement Ct — Co- Therefore, it sums up to a finite 
contribution at time-scale 
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Plugging (|2.14|) into (|2.16|) and taking expectation we obtain 



T 



(Ct) = Co - Me / dt dx X (x) d x (M(W e 0r,t))) 



T 



Co + Me dt dx X '(x) (M x {W e ix,t))), 



where we used that x is a left eigenvector of C. By evaluating as before (|2.15[) 
with S £ replaced by the true Dirac's <5-function and recalling the form (|2.13l) of the 
nonlinearity, 

1 f 1 a+ f* a f 1 a >! \( £aa(x,s) + 2K, AB ix,s) 



(Ct) « Co + eTM - dt ds J dx X 'ix) . ^ s) _ ^ s) 



where 



K{x,t) = / dyG v ix,y,t)X(p;y,y)G v (x,y,t) T . (2.17) 



The drift v — v{/3, r) is thus given by 

v = fr Hdt Cdxx'ix) ( ^ AAi fl +2 ^ Bi fl) • (2-18) 
Jo Jo \-/C B Bix,t) - 2K A Bix,t)J 

3. BOUNDEDNESS OF THE DRIFT 

The computations leading to the formula (|2.18j) for the drift have been carried 
out informally, pretending the limits we took did exist. In fact, by looking at 
(|2.18|) . it is not at all obvious that the time integral on its right hand side is finite. 
Without discussing in detail the limiting procedure performed before, in this section 
we show that such time integral is meaningful. There are two potential problems: 
the singularity of the kernel K, for t — > 0, which is the effect of the ultraviolet 
singularities mention before, and the convergence of the integral at infinity. 

We start with the analysis of the singularity around the origin. For t small 
Gix,y,t) behaves as the heat semigroup, i.e., 

G(x,y,t) =ptix- y) ^ + U(x, y, t), 

where 

^-^ 1 Qc-y+n) 2 

Ptix -y) = 2^ -n=r e 



and Hix,y,t), H y ix,y,t) are bounded as t — > + . Therefore, by (I2.17|) . the inte- 
grand on the right hand side of (|2.18l) should have the non integrable singularity 
£-3/2 f or ^ c i ose to 0. As we next show, this divergence disappears due to a can- 
celation. We compute the (dangerous) contribution to the drift coming from the 
heat semigroup part of G(x, y,t), i.e., 



with 



£ s ix,t)=[ dy [dyptix - y)} 2 S(p;y,y) = J dy g t {x - y)Y.(p\ y, y), 
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where the even function 

, \ (z + n)(z + k) (z+k) 2 +(z+ n ) 2 

behaves like (4y / 7r) _1 (5(z) as f —> + . Therefore, 

1 



l E(p;x,x)+0(t) 



so that 



s -fl f°°Ai— [ X A 'I \ \ 1 f ^aa(P;x,x) + 2Z AB (p;x,x) 
V ~ Pl Jo & 2 Jo [ ) UV* \-X BB (p;x,x)-2X AB (p;x,x) 



0(t) 



On the other hand, 



f dx '(x) ( ^aa(p;x,x) +2J: AB (p;x,x) \ 
o XX[X> \-E BB (p;x,x) - 2Y, AB (p;x i x)J 

= dx p' B [2p A (p A - 1) + <tp A p B ] + P A [2p B (p B - 1) + ^PaPb] 
Jo 

= 2 dy (p 2 A p B + p A p% - p A p B )' = 0, 
Jo 

i.e., the coefficient in front of the non integrable singularity t~ 3 / 2 vanishes. 

We now analyze the convergence of the integral (|2.18p at infinity. We notice 
that, by our assumptions on the spectrum of C, 

G(x, y, t) = 7 p'(x) x(y) + G(x, y, t), 

where the 2x2 rank one matrix r Yp~'(x)x{y) is the projection on the null space 
of C 1 while the remainder G(x,y,t) is exponentially small as t — > oo. Therefore, 
plugging this decomposition in (12.17[) we get, 

fC(x, t) = K p'(x)p'{x') T + R(x, x', t), (3.1) 

where k — j 2 J Q dy y, y)x'(y) T i while the remainder R(x,x',t) is inte- 

grable at t = oo. 

By (|2.13p . inserting (13. 1|) in (|2.18l) and integrating by parts, we deduce that the 
drift v is finite provided that 



(3.2) 



(p'bY-Wb 

As x = (Pbi ~Pa)i the above conditions reads, 



3 fdx[{p' A fp' B+ p' A {p' B f]=Q. 
Jo 



To prove (|3.2[) . we denote by I the integral on the left-hand side. Recalling that 
pc = 1 — p A — Pb and integrating by parts, we then get 

I = -3 / da; p A p B p c = dx [p A (p B p c )' + P B {p'aPcY + Pc(p'aPb)'] 

r Jo (3.3) 

= dx [(PaPbYp'c + (p a Pc)'Pb + (PbPc)'Pa]- 
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On the other hand, (|2.7|) in term of the triple (Pai Pb> Pc) reads, 

Pa = -P[{PaPb)' - (PaPcYI 
Pb = -PKPbPcY - (PbPa)% 

fc = -P[lPcPA)'-(PcpB)% 

Substituting the above relations in (|3.3[) the identity 1 = follows. 

4. Antisymmetry of the drift 

In this section we show that, as thought by Bodineau and Derrida [SJ § 6], the 
drift is antisymmetric with respect to the exchange of the species' labels. As a 
matter of fact, we present two independent arguments. The first relies on a direct 
analysis on the expression (|2.18l) , while in the second we exploit a symmetry of the 
underlying microscopic dynamics. 

Macroscopic computation. On the set of profiles p = {^f^j we introduce the 
involution O defined by 

In the language of high energy physics one may interpret O as a CP symmetry 

We then observe that if we include the dependence on the total densities r a , tb 
in the notation, the symmetry 

Op(rA,r B ) = p(rB,r A ) 
holds modulo a translation on the torus (look at (|2 - 7J> ) . 



Denoting by C the linear operator defined as C but with the matrix function B 
replaced by 

(2(Qp) A + 2{Qp) B -l 2(Qp) A \ 

V -2(0/0* -2{Qp) A -2{Qp) B + l)- [iA) 

it is easy to see that O o C = C o 0. This implies that if i/j(x,i) is solution to 
dtt/j = Cip then ip(x,t) = (Qip (-,{)) (x) is solution to d t ip = Ctp. In particular, 

G(x,y,t)=( C ! }\ G(—x, -y, t) ^ ' ' 



i o/ v ' a ' ; vi o 

where G(x,y,t) = e tc (x,y) denotes the fundamental solution associated to L. On 
the other hand, recalling (I2.6p . 



-(/':.'/■ .'/! = -(O/*: //.//! = ( ^ * j £(p; -y, -y) L () 



IC{x,t) = ( ° z) jC(-x,t) ° ' 



whence, by (|2.17p . 

,1 oy v ' y vi o 

where 

K(x,t) = / dyG y (x,y,t)T,(p;y,y)G y (x,y,t) 
Jo 
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Therefore, by (|2~T8l) 

v = pj dt dxx(x)[ ~ \ ; P \ ; 

Jo Jo y-lCAA[-X,t) - 2K.AB{-X,t) / 

= -h f°df fdx (Ox') ( X ) ( y + f~-(f ) ] (4.2) 

JO Jo y-K.BB[X,t) - 2K.AB{X,t) 

a [°° a+ f 1 A ~'f \ ( K-AA(x,t) + 2K. AB (x,t) \ 
= -h dt / dxx{x) ~ ~ ' , 

Jo Jo y-K. B B(x,t) - 2ICab(x, t)J 

where x = (O/Og, — Gp^) is the left eigenvector of C with zero eigenvalue and the 
identity \' = 0x' is straightforward. Recalling <dp(r A ,rB) = p{rB,r A ) holds mod- 
ulo a translation on the torus, by (|4.2p we conclude that f(rvi,rB) = — v(rs, rvi). 
Note indeed that by arguing as above (it is in fact easier) it holds 7(7^4,7-5) = 
l{rB, r A ). 

A microscopic symmetry. From a mathematical point of view, the ABC model 
is a continuous time Markov chain on the state space 51 = {A, B,C} Zl , where 
Zj, = {0, . . . , L — 1} is the ring of the integers modulo L. Given £ S 51, the species 
at site i is therefore ((i) € {^4, -B, C}. For /? g R the dynamical rules are specified 
by the generator \fi that acts on observables / : 51 — >• R as 

i^/(o= E c f(0[/(c" +1 )-/(o], (4.3) 

where is the configuration obtained from £ by exchanging the species at sites 

i and i + 1 and the jump rates cf are given by 

= J exp{ ^ } if (cw ' c(l + 1)} e c) ' (c ' B) ' (B ' r 4 4) 

\exp{-A} if (((i),(( i + l)) e {(4,B),(B,C),(C,i)}. 

We next prove a symmetry property of the microscopic stochastic dynamics with 
respect to a suitable involution defined on the state space 51. The antisymmetry 
of the macroscopic drift will be then shown to be a consequence of such symmetry. 
Let t be a transposition of the species' labels and denote by to <G {A, B, C} the 
image of a € {^4, B, C}. We associate to each r the involution G r : £1 — > 51 defined 

by 

e T c (*) = <(-*), 

which induces a natural transformation on observables / : 51 — > R by setting 
©r/(C) = /(@tC)' Wc claim that the microscopic dynamics satisfies, 

6 r oL-^ = l/o6 r . (4.5) 

To prove the relationship P~5)l . observe that (a, a') € {(A C), (C, B), (B, A)} if 
and only if (ra, ra') G {(A, B), (B, C), (C, A)}. Whence, by (03}, 

^(6 T C) = ^ (i+1) (C). 

Recalling (|4.3[) we deduce, 

e r oL^/(c)= X! ^ (4+1) (C)[/((e r c) M+1 )-/(e r c)]. 
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Therefore, as (6 T C) M+1 = 6 T C~ (m) '~ l , 

e T oL-^/(C)= J2 ^ (J+ i)^[/(e.r (,;+1) ^)-/(e.O] 

= E Cj(0[f(&rC" + 1 )-f(QrO] 

= L^oe r /(c). 

The relationship (|4. 5[) implies the corresponding symmetry on the statistics of the 
paths of the ABC model. In particular, if the process has a macroscopic drift v = 
v(P;r A ,r B ,r c ) then it necessarily satisfies v((3;r A ,r B ,r c ) = v(~f3;r TA ,r TB ,r TC )- 
The antisymmetry of the drift with respect to the transposition of the species' 
labels now holds provided v is an odd function of j3. Recalling that v is given by 
the expression (|2 . 1 8[) . this can be directly verified. 

Appendix A. Nonexistence of traveling waves 

In this appendix we prove the nonexistence of traveling wave solutions to the 
system (|l.ip that satisfy the constraint p A + Pb + Pc = 1. Denoting by c its speed, 
such a wave is a function of the form 

/ Y A {x-ct) \ 

p{x,t) = Y B {x~ct) 

\l-Y A (x-ct)-Y B (x-ct)J 

where Y — (Y A , Y B ) T is a noncostant one-periodic solution to 

-cY' = Y" + pN{Y)\ (A.l) 

where M( ) is defined in (|2.4p . By integration, it follows that Y solves the first 
order planar system 

Y' = F K (Y), 

where Fk(Y) = K — cY — f3JV(Y) for some constant vector K = (K A , K B ) T . Since 
divi*V(y) = —2c, the Bendixson's criterium, see, e.g., [HI § II. 4. 67], implies that for 
c 7^ the above planar system does not admit periodic (in particular one-periodic), 
non-stationary solutions. Indeed, assume by contradiction that such solution exists, 
denote by A the bounded region delimited by its orbit, and let n be the outer 
normal to A. By the Gauss-Ostrogradskij's formula, J A divFK — §q^Fk ■ n = 
since Fk -L n on d A. This yields the desired contradiction for c ^ 0. 

Appendix B. Howto guess fluctuating hydrodynamics 

We here briefly discuss, without even attempting any mathematical justification, 
how the fluctuating hydrodynamics (|2. 1[) — (|2.3[) can be inferred from the underlying 
microscopic Markovian dynamics. We refer to (5j [6j [10] for a somewhat alternative 
argument. 

Recalling the notation for the ABC dynamics introduced in Section [H the occu- 
pation numbers of the species are given by a a (i) — a G {A, B, C}, i € Til- 
From these variables we construct the empirical densities of the species as 

TT a = — ^2 a a (i)S i/L) a£{A,B,C}, 
iez L 

that we regard as a random measure on the unit torus T. 
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To infer the evolution of the empirical densities, we recall that if / : O — > R 
is an observable then, as follows from the theory of Markov chains, its expected 
infinitesimal increment is Lr /, i.e., conditionally on it holds 

(f(C(t + dt)) - f(C(t))) = L*/(c(t))d* + o(dt). 

In order to perform this computation for the empirical densities it is convenient to 
introduce smooth test functions J a : T — ► R. By setting 

(tt q , j a ) = i J a (i)v a (i), 

simple computations then show, 

iAtt q , J a ) = ~ £ [M*$r) -^a(x)]W*)[e~^o-a+i(* + l)+e^(7 a+2 (* + l)] 
iez L 



\e^a a +i(i) + e (J a+2 {i)](T a (i + 1)}, 



where the summation in the species labels a is modulo three. By approximat- 
ing discrete gradients with continuous derivatives, expanding the exponential, and 
summing by parts, we then obtain 

l/<7r«,Ja)«^ £ |4'(i) ( 7 Q+1 ( i )-/3j;(f)a C[ ( l )[a ct+1 (z + l)-a ct+2 ( l + l)]}. 

Provided that we could replace the local product a a (i)a a >(i + 1) above with the 
corresponding product of the local densities, we then identify the drift term in 
the fluctuating hydrodynamic equation (|2.1j) . While this factorization hypotheses 
has been widely used, see e.g., [6], we mention that it could be rigorously justi- 
fied by the methods of the hydrodynamical limits. More precisely, if we consider 
an initial configuration of the species that is associated to some density profile 
p(0) = (pa(0), Pb(0), Pc(0)) then the empirical densities at times 0(L 2 ) converges 
as L — > oo to the solution of the hydrodynamic equation (jl.ip with initial datum 
p(0). The proof of such a statement should be achieved by using standard tools 
in hydrodynamical limits, see e.g., [TBI [21]; however, to our knowledge, the actual 
details for the ABC model have not been carried out. We emphasize that in order 
to justify the use of the fluctuating hydrodynamics for the computation of the drift 
of phase fluctuations, we would need the validity of the factorization assumption up 
to times 0(L 3 ), certainly a most challenging issue in the context of hydrodynamical 
limits. 

To identify the fluctuations of the empirical densities, we recall that, according 
to the general theory of Markov chains, if /, g: SI — > R are two observables then, 
conditionally on £(t), 

[f(C(t + dt)) f(C(t)) l//(C(i))di] • [g(C(t + dt)) - g(at)) IA(C(£)H 

= [1/fg (C(t)) - f(at)) LMCW) - ff(C(*)) L"/(£(t))] dt + o(dt). 

The fluctuations of the empirical densities ir a are therefore characterized by the 
quadratic form, 

r a, a' {Ja > J a') — L ["(^"q ; J a) {^a' ; Joe' )1 , 

a, a' e {A,B,C}, 

— (tTq, ^q)L (tTq', Ja') — (tTq', Ja')^ {^a, J a) 
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where as before J a : T — > R are smooth test functions. By straightforward com- 
putations, for a! = a, 

T a , a {J a ,J a ) = ^ [J a (^)-j a (j;)] 2 {a a (i)[e-^a a+1 (i + l) 
iez L 

+ e&a a+2 (i + 1)] + [e^a a+1 (i) + <r a+2 (i)]a a {i + 1)} 

x {<r a (i)[<j a+1 (i + 1) + (T a+2 (i + 1)] + [a a+ t(i) + a a+2 (i)] a a (i + 1)}, 
while, for a' = a + 1, 

x {e~^a a (i)a a+1 (i + 1) + c& <7 a+1 (i)cr a (i + 1)} 
iez L 

and finally, for a' = a + 2, 

iez L 

x {e5rcr a (i)cr a+2 (i + 1) + c~ ^ (J a+2 {i)a a {i + 1)} 

w ~hh ^ J a(r) J Q +2(r){ f7 «(»K+2(* + 1 ) + ^+2W^(* + i)}. 

iez L 

In view of the factorization assumption discussed before, the fluctuations of the 
empirical densities agree with the correlation of the noise in (12.11) . Compare in 
particular (|2.2[) and (|2.3[) to the above expressions. 

Acknowledgments. We are grateful to T. Bodineau, B. Derrida, and C. Mascia 
for useful exchanges of views. 
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